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1 . 0  INTRODUCTION 

1.1  Motivation 


A  large  number  of  physical  dynamic  systems  are  nonlinear;  practi¬ 
cally  all  dynamic  systems  are  nonlinear  if  the  excitation  is  high  enough. 
The  behavior  of  most  structures,  while  subject  to  extreme  excitations, 
represents  the  design  condition.  Therefore,  accurate  knowledge  of  their 
dynamic  properties  is  crucial. 

Where  does  this  information  come  from?  For  some  systems,  the  dy¬ 
namic  properties  can  be  computed  theoretically.  In  some  limited  cases, 
specialized  subtests  can  be  conducted  and  the  results  combined  to  ap¬ 
proximate  the  total  system  behavior.  More  often,  however,  the  complexi¬ 
ty  of  the  nonlinear  interactions  make  it  necessary  to  perform  tests  on 
the  complete  system  to  determine  its  behavior.  If  the  analyst  is  ex¬ 
tremely  fortunate,  there  may  even  be  full-scale  measurements  of  excita- 
tion(s)/system  response(s)  directly  representing  the  extreme  loading 
(design)  case.  A  more  common  situation  is  to  have  one  set  of  measure¬ 
ments  representing  the  excitation/responses  for  a  level  of  excitation 
variance  lower  than  the  design  level. 

Even  when  an  accurate  data  set  is  available,  the  analyst  is  only 
slightly  better  off  than  before  because  nonlinear  system  theory  is  not 
well  developed.  A  "single  useful  method"  for  identifying  and  quantify¬ 
ing  nonlinear  system  properties  does  not  exist.  While  some  techniques 
do  exist  and  could  possibly  be  applied  by  an  experienced  user,  the 
results  are  limited,  expensive,  and  difficult  to  interpret. 

The  U.S.  Navy  often  faces  this  scenario  while  designing  new  proto¬ 
type  systems  for  which  there  is  no  historical  basis  for  design.  The 
development  of  a  general  nonlinear  system  identification  procedure  has 
been  underway  at  the  Naval  Civil  Engineering  Laboratory  (NCEL)  since 
1980,  motivated  primarily  by  the  need  for  an  improved  understanding  of 
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vessel  moorings  and  dynamic  wave  loads  (Ref  1).  Development  of  the 
basic  frequency-domain  mathematical  framework  used  as  the  foundation  of 
this  new  method  was  partially  sponsored  by  NCEI.  and  has  been  recently 
published  (Ref  2). 

This  report  summarizes  the  development  of  a  new  approach  to  nonlin¬ 
ear  system  identification  developed  through  NCEI.  that  has  been  success¬ 
fully  applied  to  a  variety  of  simulated  nonlinear  systems. 

1.2  State  of  the  Art  in  Higher  Order  System  Identification 

General  nonlinear  systems  exhibit  a  more  complex  behavior  than 
linear  systems.  Accordingly,  much  of  the  research  in  nonlinear  systems 
theory  has  been  directed  toward  defining  a  universally  applicable 
mathematical  model  capable  of  simultaneously  converting  one  arbitrary, 
stationary,  input  variable  into  the  many  terms  that  comprise  the  non¬ 
linear  output.  Such  a  theory  has  not  yet  been  discovered. 

The  nonlinear  formulation  with  the  most  universal  application  is 
the  so-called  Volterra  series  (Ref  3).  The  Volte rra  series  assumes  that 
the  nonlinear  response  can  be  expanded  as  a  "power  series  with  memory," 
with  an  N- integral  kernel  associated  with  each  nth  order  component. 

This  approach  leads  to  higher  order  spectra  such  as  bispectra  and 
tri spectra  that  are  based  on  higher  order  cumulants  of  the  input  time 
history.  Some  of  the  disadvantages  of  Volterra  series  are: 

(a)  Unreasonable  computer  processing  requirements  for  kernels  of 
order  four  or  more. 

(b)  Sensitivity  to  the  distribution  of  the  input  signal  and 
noise . 

(c)  Difficulty  in  interpreting  system  behavior  by  inspection  of 
third  and  higher  order  kernels. 

(d)  Large  random  errors  that  require  excessive  (often  unobtainable) 
amounts  of  data  for  acceptable  confidence  limits. 
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Thus,  in  practice,  the  Volterra  series  is  typically  applied  only  when 
the  highest  order  response  is  assumed  to  be  second-order  or  third-order 
It  is  noted  that  the  Volterra  kernels  are  invariant  to  input  excitation 
levels  only  when  the  assumed  and  physical  highest  order  responses  are 
the  same. 

Reference  3  provides  a  very  complete  overview  of  the  Volterra  and 
related  Wiener  nonlinear  system  theories.  Reference  4  is  a  similarly 
complete  review  of  the  theory,  computation,  and  applications  of  the  bi¬ 
spectrum.  Autoregressive  (AR),  Moving  Average  (MA) ,  or  combined  (ARMA) 
models  have  also  been  used  to  describe  higher  order  system  responses 
(Ref  5).  However,  these  models  do  not  provide  system  parameters  appro¬ 
priate  for  structural  problems.  Reference  6  summarizes  much  of  the 
latest  research  activity  and  applications  in  higher  order  spectral 
analys is . 

1.3  Problem  Description 

As  stated  in  Section  1.1,  NCEL  has  had  a  continual  need  for  nonlin 
ear  systems  identification  techniques.  Some  specific  examples  include: 

(a)  Wave- induced  forces  on  slender  members  (the  Morison  equation) 

(b)  Nonlinear  stiffness:  catenary  (mooring  lines),  fenders 
(pierside) . 

(c)  Second-order  wave  drift  forces  on  ships  (proportional  to  the 
wave  envelope). 

(d)  Extreme  vessel  motions. 

(e)  Buoy  motions  (nonlinear  buoyancy  for  floating  versus 
submerged) . 

(f)  Time-varying  stiffness:  fishtailing  in  single  point  moorings 
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(g)  Time-varying  mass  and  stiffness:  bottom  interaction  of 
mooring  lines. 

For  a  nonlinear  system  like  a  moored  ship,  the  above-mentioned 
matters  may  occur  simultaneously.  In  this  case,  if  a  full-scale  data 
set  is  available,  it  typically  consists  of  measurements  of  the  inputs 
consisting  of  the  incident  wave  and  wind  vectors,  and  the  system  res¬ 
ponses  of  interest,  usually  consisting  of  the  six  vessel  motions  and 
selected  tensions  in  the  mooring  lines  (Ref  7).  For  slender  member  wave 
force  data  sets,  the  input  might  consist  of  either  the  incident  waves  or 
wave  velocities  while  the  output  is  the  horizontal  component  of  force  in 
line  with  the  direction  of  wave  advance.  For  the  finite  duration  of  the 
data  set,  it  is  usually  assumed  that  the  input  data  are  stationary  with 
constant  variance. 

In  all  cases,  the  purpose  for  collecting  data  is  to  allow  for  vali¬ 
dation  of  a  postulated  engineering  model  of  the  system.  Usually,  this 
model  consists  of  component  terms  such  as:  inertial  forces;  linear  and 
viscous  damping  forces;  restoring  forces  that  represent  the  equation  of 
motion  of  the  system;  and  component  terms  representing  the  excitation 
forces  (for  example,  first-  and  second-order  hydrodynamic  loads).  The 
validation  seeks  to  identify  the  nature  of  the  physical  process  (qualita¬ 
tive  properties)  and  the  associated  numerical  parameters  (quantitative 
properties) . 

1.4  Overview  of  Present  Study 

The  objective  of  the  present  study  was  to  develop  a  general 
nonlinear  system  identification  methodology  to  determine  from  measured 
input  and  output  stochastic  data  the  validity  of  each  term  in  a  proposed 
nonlinear  differential  equation  of  motion.  This  report  describes  some 
of  the  results  obtained  by  a  new  practical  frequency-domain  approach 
that  has  been  under  investigation  at  NCEL  since  1980  (Ref  2  and  8), 
where  established  multiple-input/single-output  linear  procedures  are 
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applied  to  reverse  dynamic  systems.  Material  in  this  report  focuses  on 
single  deeree-of -freedom  (SDOF)  nonlinear  systems  excited  by  stationary 
random  data;  this  excitation  can  be  Gaussian  or  non-Gaussian  with 
arbitrary  spectral  shape.  The  new  approach  is  theoretically  outlined 
here  and  applied  to  a  variety  of  nonlinear  systems. 

It  is  noted  that  this  approach  was  also  used  independently  in  Ref¬ 
erences  9  and  10  based  on  ideas  in  Reference  8.  This  report  greatly 
extends  the  understanding  and  application  of  the  material  contained  in 
References  8  to  10.  It  is  demonstrated  here  that  this  new  NCEL-developed 
method  identifies  desired  nonlinear  system  properties  independent  of  the 
probability  or  spectral  nature  of  the  measured  input  and  output  data, 
excitations  levels,  or  typical  extraneous  noise  levels.  Simulated  SDOF 
nonlinear  systems  include  Duffing,  Van  der  Pol,  Mathieu,  and  Dead-Band 
systems . 
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2.0  SINGLE  DEGREE -OF -FREEDOM  (SDOF)  NONI.INF.AR  SYSTEMS 


2.1  Systran  Descriptions 

A  broad  class  of  nonlinear  SDOF  dynamic  systems  is  described  by  the 
constant  coefficient  differential  equation: 

mu  +  cu  +  kn  +  p(  u  ,  u ,  t. )  =  F(  t )  (2.1) 

For  the  mechanical  systems  addressed  in  this  study,  m  is  the  system 
mass,  k  is  the  linear  elastic  stiffness  constant,  c  is  the  linear  viscous 
damping,  and  p(u,u,t)  represents  a  general  nonlinear  damp i rig- i es tor i ng 
force  function  component.  The  external  applied  force  is  F(t)  and  the 
displacement  response  is  u(t)  or  simply  n .  From  a  mechanical  schematic 
viewpoint.,  as  described  in  Figure  1,  the  above  differential  equation 
represents  dynamic  equilibrium  following  the  relationship: 

F  ( u , u , t )  =  F( t )  -  mu  (2.2) 

r 

where  the  total  restoring  force  is 

F  ( i.i ,  u  ,  t )  =  cu  +  ku  +  p(u,u,t)  (2.3) 

The  basic  linear  case  occurs  when  p(u,u,t.)  is  zero. 

A  variety  of  nonlinear  SDOF  syste  ,s  are  considered  in  this  paper. 
They  were  selected  to  represent  different  nonlinear  behavior  as  described 
by  the  following  nonlinear  force  components: 


(a)  Duffing  system  (Ref  11): 

3 

p( n,u, t)  =  k^u 

(b)  Van  der  Pol  system  (Ref  11): 

p(",u,t)  =  -3  j[u3| 


(2.4) 


(2.5) 
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(c)  Combined  Duffing  -  Van  der  Pol  system: 


p(u,u,t)  =  jk  + 


'3  Id 


(2.6) 


(d)  Mathi eu  system  (Ref  11): 


p(n,u,t)  =  k  cos(2n  f  +  $  )  u 
m  mm 


(e)  Dead-Band  system: 

-  ku , 


(2.7) 


p(u,u, t) 


’{ 


for  u  <u 


(2.8) 


ku  sgnfu)  for  lul  >  u 
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(f)  Extreme  Duffing  -  Van  der  Pol  system: 

Case  (c)  above,  with  linear  coefficients  c 


=  k  =  0. 


All  of  the  above  types  of  nonlinearity,  except  for  the  Mathien  type, 

describe  systems  whose  characteristics  are  nonlinear  with  amplitude. 

The  Mathieu  "nonlinearity"  in  actuality  describes  a  linear  system  with  a 

periodic,  time-dependent  stiffness.  The  nonlinearity  for  the  Dead-Band 

system  is  defined  so  that  the  effective  system  stiffness  [ku  +  p(u,u,t)] 

is  zero  for  I u I <u  and  k[u-u  sgn(u)l  for  I u I >u  . 

8  8  8 

2.2  Simulation  of  Measured  Dynamic  Response 


When  a  nonlinear  dynamic  system  model  is  subjected  to  a  general 
time-dependent  excitation  F(t),  the  dynamic  response  history,  u(t),  is 
normally  determined  by  numerical  integration.  For  the  present  study, 
the  Adams  -  Bash  for th  pred ict m -corrector  method  (Ref  12)  is  employed. 

The  excitation  force  F(t)  used  in  all  cases  is  broad  band  random 
data  following  a  Gaussian  probability  distribution.  The  random  force 
history  is  formed  using  a  random  number  generator  provided  in  the  MATI.AB 
386  computer  code  (Ref  13).  The  discretization  time  At  is  0.02  seconds, 
which  defines  the  Nyquist  cut-off  frequency  at  28  Hz.  A  low-pass  filter 
is  applied  to  the  raw  excitation  history  to  form  an  applied  force  history 
with  a  nearly  uniform  spectral  content  below  12.3  Hz.  A  record  length 
of  8,192  samples  was  used  in  all  simulations. 


Dynamic  properties  for  all  simulated  systems  were  limited  to  the 
lower  frequency  range  (below  10  Hz).  This  choice  was  made  to  assure 
that  the  numerical  simulations  using  a  0.02-second  time  increment  would 
bo  relatively  error  free. 

For  situations  considered  in  this  paper  where  extraneous  noise 
effects  are  studied,  a  nroad  band  Gaussian  random  history,  uncorrelated 
with  the  excitation,  is  generated  and  added  to  the  calculated  response 
u(t).  This  represents  a  measured  displacement  corrupted  by  noise.  The 
applied  force  is  always  assumed  to  be  noise  free.  Numerical  values  for 
the  noiseless  linear  and  nonlinear  systems  in  this  paper  are  listed  in 
Table  1  in  Section  5.1. 
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3.0  PRELIMINARY  ANALYSIS  OF  MEASURED  DATA 


3.1  Objectives  and  Existing  Practice 

The  objectives  of  preliminary  data  analysis  for  a  suspecte.d  nonlinear 
sys  t  em  are : 

(a)  Recognition  that  nonlinearity  is  present  (or  absent). 

(b)  Identification  of  the  general  class  of  nonlinearity  (to 
minimize  analysis). 

(c)  Estimation  of  some  key  system  parameters  (in  some  cases). 

A  variety  of  data  descriptions  may  be  used  to  accomplish  these  ob¬ 
jectives  consisting  of  excitation  and  response  probability  densities 
(histograms),  excitation  and  response  autospectra,  and  constructed 
force/displacement  plots. 

For  the  case  of  a  linear  system  subjected  to  Gaussian  random 
excitation,  both  the  excitation  and  response  probability  densities  are 
theoretically  Gaussian.  However,  if  the  system  is  nonlinear,  the  res¬ 
ponse  probability  will  deviate  from  the  Gaussian  form,  often  in  a  manner 
suggesting  the  nature  of  the  nonlinearity  (Ref  2). 

A  linear  SDOF  system  that  is  excited  by  broad  band  random  data 
typically  exhibits  response  in  a  narrow  band  in  the  vicinity  of  the 
resonant  frequency.  When  nonlinearity  is  present,  the  response  of  the 
system  may  contain  higher  harmonics  relative  to  the  resonant  frequency 
band.  This  effect  may  be  noted  in  the  response  autospectrum .  In  addi¬ 
tion,  for  the  case  of  a  Math  leu-type  system,  it  will  be  shown  that  the 
response  autospeetrum  provides  an  estimate  of  the  periodic  stiffness 
f  requency . 
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3.2  Restoring  Forrc/Displacement  Relationships 

The  construction  of  force/displacement  relationships  offers  an 
informative  means  of  uncovering  the  general  nature  of  a  subject  dynamic 
system.  Given  the  measured  (or  calculated)  acceleration  response  u(t) 
and  applied  force  F(t),  the  time  history  of  the  SDOF  system  restoring 
force  Fr(u,u,t)  can  be  constructed  provided  the  system  mass  m  is  known. 
This  calculation  follows  Equation  2.2.  It  will  be  shown  in  the  following 
discussion  that  the  force/displacement  patterns  associated  with  response 
to  random  and  harmonic  excitation  have  similar  characteristics  whether 
the  system  is  linear  or  nonlinear.  These  patterns  clearly  reveal  the 
general  nature  of  the  system  nonlinearity.  Patterns  for  linear  systems 
are  presented  first,  followed  by  each  of  the  nonlinear  systems  described 
in  Section  2.1. 

Consider  a  linear  SDOF  system  with  viscous  damping  undergoing  simple 
harmonic  motion  (SHM)  where: 

u(t)=Usin(2irft)  (3.1) 

o 

The  corresponding  restoring  force  history  (see  Equation  2.3  with 
p(u,u,t)=0)  is: 

F  (u,u,t)  =  k  U  sin(2nf  t)  +  (2irf  )  c  U  cos(2nf  t)  (3.2) 

r  o  o  o 

This  pair  of  functions,  when  plotted  with  u(t)  as  the  abscissa  and 
Fr(u,u,t)  as  the.  ordinate,  will  trace  a  closed  oval  as  illustrated  in 
Figure  2a.  When  the  same  dynamic  system  is  subjected  to  broad  band 
random  excitation,  a  similar  geometric  pattern  having  the  same  oval 
envelope  is  obtained  as  illustrated  in  Figure  2b.  This  similarity 
between  harmonic  and  random  excitation  based  force/displacement  plots 
occurs  since  the  response  of  a  linear  SDOF  system  is  narrow  band;  the 
general  appearance  of  the  displacement  response  is  that  of  a  fluctuating 
sinusoid  with  frequency  close  to  the  SDOF  natural  frequency.  The  slope 
of  the  mean  value  of  F^  at  any  given  value  of  u  corresponds  to  the  linear 
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stiffness  k.  In  addition,  the  width  of  the  oval  at  u=0  Indicates  the 
value  of  (2TTf^)cU.  This  pattern  description  applies  to  both  simple 
harmonic  motion  and  random  motion. 

For  the  nonlinear  systems  under  investigation,  the  patterns 
associated  with  restoring  force  versus  displacement  exhibit  a  property 
similar  to  that  noted  for  linear  systems.  In  particular,  the  closed 
curve  traced  by  imposing  a  simple  harmonic  displacement  history  on  the 
appropriate  restoring  force  relationship  (combining  Equations  2.4  to  2.8 
with  2.3)  produces  a  characteristic  nonlinear  "signature."  The  nonlinear 
force/displacement  pattern,  due  to  broad  band  random  excitation,  traces 
a  pattern  whose  envelope  is  similar  to  that  associated  with  the  simple 
harmonic  displacement  based  pattern.  Force/displacement  patterns  corres¬ 
ponding  to  the  example  nonlinear  systems  subjected  to  both  types  of  ex¬ 
citation  (with  system  properties  listed  in  Table  1)  are  described  below: 

(a)  The  Duffing-type  restoring  force  versus  displacement  traces  an 
envelope  with  a  mean  slope  following  the  stiffness  nonlinearity  as 
depicted  in  Figure  3.  The  hardening  spring  behavior  expected  from  the 
Duffing  (cubic)  restoring  force  component  is  clearly  evident. 

(b)  The  Van  der  Pol-type  restoring  force  versus  displacement 
traces  an  envelope  with  a  mean  slope  following  a  constant  linear  stiff¬ 
ness.  The  bulbous  envelope  at  the  larger  displacement  levels,  shown  in 
Figure  4,  is  indicative  of  the  Van  der  Pol  damping  nonlinearity. 

(c)  The  combined  Duffing-Van  der  Pol-type  restoring  force,  shown 
in  Figure  5,  traces  an  envelope  which  exhibits  the  stiffness  and  damping 
non  1  inear  it ies  in  a  manner  consistent  with  each  of  the  behaviors  described 
above . 

(d)  The  Mathieu-type  restoring  force  versus  displacement  pattern, 
shown  in  Figure  6,  traces  a  fluctuating  envelope  which  alternates  between 
a  low  and  high  value  of  "linear"  stiffness. 
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(e)  The  Dead-Band  (or  Gapping) -type  system  exhibits  a  clear  envelope 
of  restoring  force  versus  displacement,  shown  in  Figure  7,  which  enables 
the  investigator  to  estimate  the  value  of  the  dead  band  (in  this  example, 
between  +0.0025).  In  a  later  section  it  will  be  noted  that  this 
preliminary  dead-hand  estimate  is  required  for  identification  of  this 
type  of  nonlinear  system  using  the  spectral  method. 

3.3  Further  Preliminary  Analysis  for  Mathieu-Type  Systems 

If  a  SDOF  system  is  suspected  to  be  of  the  Mathieu-type  on  the  basis 

of  the  restoring  force/displacement  plot,  additional  preliminary  analysis 

may  be  performed  to  estimate  the  frequency  of  the  fluctuating  stiffness 

(f  ).  It  can  be  demonstrated  analytically  that  the  response  of  such  a 

system  to  broad  band  random  excitation  contains  significant  components 

at  the  "mean"  natural  frequency  f  =  ( l/2ir)/k/m,  and  components  at  the 

absolute  sum  and  difference  frequencies  I f  +f  |  and  If  -f  |.  This 

1  n  m‘  1  n  m 1 

property  is  clearly  illustrated  in  the  displacement  response  auto¬ 
spectrum  shown  in  Figure  8  for  a  system  with  f  =  3  Hz  and  f  =  1  Hz. 

n  m 

Some  caution  must  be  taken  in  the  interpretation  of  displacement 

response  autospectra,  especially  for  cases  in  which  f  >  f  .  For 

m  n 

example,  if  f  =  3  Hz  and  f^  =  5  Hz,  then  significant  frequency 
components  will  occur  near  2,  3,  and  8  Hz. 
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4.0  NONLINEAR  SYSTEM  IDENTIFICATION  TECHNIQUE 


4. 1  Reverse  Dynamic  Systems 

The  spectral  analysis  technique  developed  in  this  paper  for  identifi¬ 
cation  of  SDOF  nonlinear  systems  makes  use  of  a  "reverse"  dynamic  system 
viewpoint  (i.e.,  a  reversal  of  the  input/output  roles  of  excitation  and 
response).  The  advantage  of  this  approach  is  noted  by  examining  a  non¬ 
linear  system  the  usual  (forward)  way.  The  excitation/response  diagram 
is  shown  in  Figure  9  for  a  general  SDOF  nonlinear  system  described  by 
Equation  2.1.  The  linear  frequency  response  function  defined  by: 

H  (  f )  =  [k  -  (2trf  )2m  +  j  (2irf)c]_1  (4.1) 

relates  an  "effective"  excitation  force  input  Fg(t)  versus  displacement 
response  output  u(t),  where  the  effective  force  is  defined  as: 

Fg(t)  =  F(t)  -  p(u,u,t)  (4.2) 

with  p(u,u,t)  as  a  feedback  term  and  F(t)  as  the  applied  external  excita¬ 
tion  force.  Due  to  the  presence  of  the  feedback  term,  identification  of 
the  nonlinear  system  requires  use  of  an  (awkward)  iterative  approach 
when  the  above-noted  excitation/response  viewpoint  is  taken. 

Consider  now  the  alternative  viewpoint  in  which  the  displacement 
response  u(t)  is  considered  as  the  input  in  Equation  2.1,  and  the  excita¬ 
tion  force  F(t)  required  to  produce  that  response  is  calculated  as  the 
output.  This  reverse  dynamic  system  is  shown  in  Figure  10.  The  non¬ 
linear  SDOF  system  viewed  in  this  manner  is  now  a  multiple-input/single¬ 
output  (MI/SO)  system  where  the  number  of  multiple  inputs  depends  on  the 
nature  of  the  nonlinearity  term  p(u,u,t). 

All  of  the  nonlinear  system  types  treated  in  this  paper  can  be 
described  by  the  generic  reverse  dynamic  three- input/s  ingle-output 
linear  model  in  Figure  11.  The  input/output  path  A ^ ( f )  is  simply  the 
reciprocal  of  the  linear  frequency  response  function  H(f),  namely, 


(4.3) 


Aj(f)  =  k  -  (2irf)2m  +  j(2irf)c 


2 

+  j  2C 


w  i  t  h 


f  =  kz  ^  k/m 

n  2  IT 


(4.4) 


and 


e 


c 

4tt  f  m 
n 


(4.5) 


The  nonlinear  paths  are  modeled  using  the  new  frequency  domain 
technique,  described  in  Reference  2,  that  separates  each  nonlinearity 
into  separate  amplitude  and  frequency  functions.  This  technique  is 
fully  described  in  the  next  section. 

For  each  system  type  in  Section  2.1,  the  nonlinear  path  functions 
are  described  below,  where  g^(u)  represents  a  zero-memory  (amplitude 
domain)  nonlinear  function  in  Figure  11  and  A^(f)  represents  the 
associated  linear  frequency  response  function: 


(a)  Duffing 

g2(u)=u3  A2(f)=k3  g3(u)=0  A3(f)=0  (4.6) 

(b)  Van  der  Pol 

a  c3 

g2(u)=ti  A2(f)=j(2nf)  y  g3(")=0  A3(f)=0  (4.7) 

(c)  Combined  Duffing-Van  der  Pol 

a  ca 

g2(u)=u“  A2(f)=k3  +  (2lTf)  3  83(u)=0  A3(f)=0  (4.8) 

(d)  Mat'll  ieu 

g  ( u )=u  cos(2irf  t.)  A  (f)=  k  cos(0  )  (4.9) 

2  m  2mm 


(e)  Dead-Band 


82(u)  = 
83(«)=0 
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u  sgn(u) 

O 


for  Ini <u 

g 

for  1 11 1  >u 

g 


A3(f)=0 


A2(f)= 


(4. 10) 


By  defining  the  subject  nonlinear  dynamic  systems  in  terms  of  Equa¬ 
tions  4.3  to  4.10,  the  unknowns  consist  of  the.  A  ^  ( f )  frequency  response 
functions.  The  known  data  are  u(t),  F(t),  and  the  constructed  histories 
defined  by  g^ (u) . 

Using  the  new  spectral  analysis  technique  and  reverse  dynamics  con¬ 
verts  the  original  single-input/single-output  (SI/SO)  nonlinear  problem 
into  an  equivalent  multiple-input/single-output  (M1/S0)  linear  problem. 
Conversion  of  the  original  ,  di.ff  icul  t-to-analyze  nonl  inear  problem  into 
an  easy-to-ana lyze  linear  problem  is  a  significant  feature  of  this  tech¬ 
nique.  The  general  random  data  analysis  procedures  for  identification 
of  linear  frequency  response  functions  A ^ ( f )  developed  in  References  14 
and  15  are  presented  next.  Note  that  the  g.  functions  can  represent  a 
wide  range  of  physical  non  1  inear it ies  beyond  the  sample  cases  examined 
in  this  report. 


4.2  MT/SO  Analysis  of  Reverse  Dynamic  Systems  with  Zero-Memory 
Non  1 inea r i t ies 


The  usual  SDOF  linear  system  satisfies  the  linear  differential 
equation  of  motion: 


mu(t)  +  cu(t)  +  ku(t)  =  F(t) 


(4.11) 


wher° 


f(t)  =  force  input 
u(t)  =  displacement  output 
m  =  mass  of  the  system 
c  =  linear  damping  coefficient 
k  =  linear  stiffness  coefficient 
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The  physical  parameters  m,  c,  and  k  are  assumed  to  be  constants. 

For  linear  systems,  a  Gaussian  force  input  produces  a  Gaussian  displace¬ 
ment  output. 

A  SDOF  system  becomes  nonlinear  when  there  are  nonlinear  stiffness 
terms  and/or  nonlinear  damping  terms.  With  nonlinear  stiffness  terms  as 
an  example,  a  typical  nonlinear  differential  equation  of  motion  is: 

mu(t)  +  cu(t)  +  ku(t)  +  k^u^Ct)  +  k^u3(t)  =  F(t)  (4.12) 

2  3 

where  k^  and  k^  are  nonlinear  stiffness  coefficients  for  u  and  u  , 
respectively.  Here,  a  Gaussian  force  input  produces  a  non-Gaussian 
displacement  output. 

To  apply  the  nonlinear  spectral  analysis  techniques  developed  in 
Reference  2,  the  roles  of  F(t)  and  u(t)  are  reversed  as  rreviously  dis¬ 
cussed.  Equation  4.11  then  becomes  for  the  linear  case: 

y  C  t )  =  mx(t)  +  cx(t)  +  kx(t)  (4.13) 

where  x(t)  =  model  input  =  physical  output  =  displacement  u(t) 
y(t)  =  model  output  =  physical  input  =  force  F(t) 

For  cases  of  nonlinear  stiffness  terms  governed  by  Equation  4.12, 
when  the  roles  of  F(t)  and  u(t)  are  reversed,  Equation  4.12  becomes: 

y(t)  =  mx( t )  +  cx(t)  +  kx(t)  +  k2x2(t)  +  k3x3(t)  (4.14) 

It  is  not  necessary  here  to  assume  that  the  input  x(t)  is  Gaussian; 
the  technique  is  applicable  to  arbitrary  data. 

Fourier  transforms  of  both  sides  of  Equation  4.13  yield  the  linear 
frequency  domain  formula: 

Y ( f )  =  A1(f)X1(f)  (4.15) 
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where  Y(f)  =  -J.  [y(t)} 


(4.16) 


Xj(f)  =  ;}  [x(t)}  (4.17) 

A  (f)  =  k  -  (2irf)2m  +  j(2irf)c  (4.18) 

Note  that  A^(f)  Is  the  reciprocal  of  the  usual  linear  frequency 
response  function  H(f)  for  Equation  4.11  as  given  in  Equation  4.1. 

For  the  nonlinear  equation  with  constant  coefficients  case  of 
Equation  4.14  (zero-memory  nonlinearities),  Fourier  transforms  of  both 
sides  yield  the  nonlinear  frequency  domain  formula: 

Y ( f )  =  A1(f)X1(f)  +  k2X2(f)  +  k3X3(f)  (4.19) 

where  Y(f),  X  (f),  and  A .(f)  are  given  by  Equations  4.16  through  4.18 

^  ^  2  3 

and  where  X^(f)  and  X^(f)  are  the  Fourier  transforms  of  x  (t)  and  x  (t) 

denoted  by: 

X?(f)  =  $  (x2(t)j  (4.20) 

X3(f)  =  ;}  (x3(t)  j  (4.21) 

Equation  4.19  can  be  displayed  as  a  SI/SO  nonlinear  model  as  shown 
in  Figure  12. 

Figure  12  can  be  reconfigured  into  an  equivalent  three-input/single- 
output  linear  model  with  correlated  inputs  X^(f),  X2(f),  and  X^(f)  as 
shown  in  Figure  13.  They  produce  correlated  outputs  denoted  by  Y^(f), 
Y2(f),  and  Y3(f). 

Practical  procedures  to  identify  A^(f),  k2>  and  from  measured 
random  data  of  x(t)  and  y(t)  are  developed  in  References  14  and  15. 

These  procedures  show  how  to  replace  the  given  set  of  correlated  inputs 
with  a  new  set  of  uncorrelated  inputs  that  are  calculated  using  condi¬ 
tioned  spectral  density  functions.  There  is  no  requirement  that  the 
input  data  be  Gaussian. 
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4.3  MI/SO  Analysis  of  Reverse  Dynamic  Systems  with  Finite-Memory 
Non linearities 

More  extensive  types  of  nonlinear  stiffness  properties  can  occur  in 
physical  problems  when  the  nonlinear  terms  have  finite  memory  and  the 
system  operators  and  become  arbitrary  functions  of  frequency 
representing  finite-memory  convolutions  in  the  time  domain.  In  these 
cases,  Equation  4.19  takes  the  more  general  form: 

Y(f)  =  AjCflXjCf)  +  A2( f )X2(f)  +  A3(f)X3(f)  +  N(f)  (4.22) 

where  N(f)  represents  extraneous  uncorrelated  output  noise.  Figure  13 
then  becomes  Figure  14. 

For  arbitrary  input  data  x(t),  by  appropriate  processing,  Figure  14 
can  be  replaced  by  the  equivalent  model  of  Figure  15  with  new  uncorrelated 
inputs  X ^ ( f ) ,  ^(f),  and  2j(f)  passing  through  new  linear  systems 

LjCf),  L2(f),  and  L3(f). 

The  uncorrelated  inputs  in  Figure  15  can  be  denoted  by  the  simpler 

notation  Z^(f)  as  shown  in  the  figure.  These  inputs  produce  uncorrelated 

outputs  denoted  by  Y  (f),  Y,(f),  and  Y  (f)  as  also  shown. 

9  D  G 

Relations  for  Figures  14  and  15  are  as  follows  using  the  notation 
in  Reference  15  for  various  conditioned  spectral  density  functions: 


h  j(f)=  Gly(f)/Gn(f)  (  4 . 23) 

L2(f)=  G2yl(f)/G22-l(f)  (4'24) 

h3(f)=  G3y.2!(f)/G33.2!(n  (4’25) 

Also, 

A3(f)=  L3( f )  (4.26) 

A2(f)=  L2(f)  -  [G23. 1(f)/G22. 1(f)]A3(f)  (4.27) 


Aj(  f  )=  Lj(f)  -  [G12(f)/6n(f)]A2(f)  -  [G13(f)/Gn(f)]A3(f)  (4.28) 
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The  system  I.j(f)  is  the  optimum  linear  system  to  obtain  y(t)  from 
x(t).  In  general,  A  ^  (  f )  f  I.^f). 

4.4  Goodness -of-F it  via  Coherence  Functions 


It  is  standard  practice  in  the  analysis  of  linear  MT/SO  systems  to 
calculate  various  coherence  functions  to  quantify  the  goodness-of - f i t  of 
the  chosen  stochastic  model  and  the  numerical  results.  This  is  a  vitally 
important  step  because  it  provides  an  absolute  measure  of  how  well  each 
term  of  the  postulated  model  fits  the  data  over  specific  regions  (i.e., 
each  frequency).  This  level  of  information  is  not  readily  available  in 
time-domain  system  identification  techniques. 

Tn  general,  care  must  be  taken  when  coherence  functions  are  calculated 
for  correlated  inputs.  In  these  cases,  the  recommended  practice  is  to 
uncorrelate  the  inputs  via  conditioning,  then  compute  ordinary  coherences 
which  can  be  inspected  individually  and  summed  algebraically  to  get  the 
goodness-of-f i t  for  the  entire  model.  This  was  the  reason  for  reconfigur¬ 
ing  the  problem  into  the  form  shown  in  Figure  15  with  uncorrelated  inputs 

Z.(f)  and  correspondingly  uncorrelated  outputs  Y  (f),  Y  ( f ) ,  and  Y  (f). 
j  D  C 

Further  details  regarding  proper  calculation  and  interpretation  of  co¬ 
herence  functions  follows. 


The  ordinary  coherence  function  21(f)  states  the  proportion  of  the 
output  spectrum  G  (f)  that  is  due  to  X^(f)  passing  through  the  linear 
system  I.^(f)  that  gives  all  the  ways  X ^  ( f )  can  go  to  the  total  output 
Y(f).  This  quantity  is  the  usual  1 inear  coherence  function  between  the 
input  Xj(f)  and  the  output  Y(f)  as  defined  by: 


=  *Vf)  =  |Gly(f,|2/!cxx<f)  Vf)1 


(4.29) 


2 

The  partial  coherence  function  3f  ^  ( f )  states  the  proportion  of 
the  conditioned  output  spectrum  G  .(f)  that  is  due.  to  the  conditioned 

yy.J 

input  X^  j(f)  passing  through  only  the  system  to  get  the.  conditioned 
output  Yy  ^(f)  defined  by: 


*2y  l,r>  -  IV,<f)|2/lG22.,(f)  Gyy-l(r)l 


(4.30) 
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The  nonlinear  coherence  function  q  (f)  from  Reference  2  states 

-  f 


the  proportion  of  the  total  output  spectrum  G  (f),  not  ^yy  j ( f ) »  that 
is  due  to  the  conditioned  input  ^(f)  passing  through  only  the  1,^ 
system  to  go  to  the  total  output  Y(f),  not  Y^  ^ ( f ) .  This  is  defined  by: 

2  /  r  \  I 


xYw 


(f)  = 


lG2yl(f)|  /lC22-i(f>  Gyy(,)l 


(4.31) 


But 


yy 


•  i<f)  =  [i  -  »lycnl  Gw(f) 


yy 


(4.32) 


Hence 


i...  <n  -  I'  ■  4<f)) 


xyt 


2  y  1 


(4.33) 


The  next  partial  coherence  function  2!^^  states  the  proportion 

of  the  conditioned  output  spectrum  G  2!^^  that  is  due  to  the 

conditioned  input  2|(f)  passing  through  the  system  to  get  to  the 

conditioned  output  Y  (f),  defined  by: 

y  ■  z . 


,„(0  -  ISy.2.(r,|2/[G„.,,(f)  G„v., ,(f)] 


3y  •  2 ! 


33-2! 


yy  ■  2 ! 


(4.34) 


The  nonlinear  coherence  function  q  (f)  in  Reference  2  states  the 

’xy 

c 

proportion  of  the  total  output  spectrum  G  (f),  not  G  2!^^’  is 

due  to  the  conditioned  input  X^  ^ ^ )  Passing  through  only  the  system 

to  get  to  the  total  output  Y(f),  not  Y  ^(f).  This  is  defined  by: 


4  <r)  =  !fi3y.2!(f)|2dG33.2,(n  Gyy(f)] 


(4.35) 
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Rut 


r,  (f)  =  i  -  y;  .(f) |  g  ,(f) 

yy • 2 !  1  2  y  -  1  1  yy • 1 


*2y  l(f)l  I'  -  *Vn|  V0 


( A . 36) 


Hence  , 


q2  (n  =  fi  -  y2  ,(f)|  fi  -  y2  m|  y2  olm 

'xyr  1  2y -  1  11  ly  1  By -2! 


(4.37) 


With  these  component  cohprpncp  f  unctions  defined,  t  Iip  summed 
coherence  function  for  the  pntirp  model  r.an  1>p  constructed.  Using 
Rpfproncp  2,  tin's  is  defined  as  tlie  multiple  coherence  function  and  is 
given  for  this  three- input  nonlinear  model  as: 


y2  (f)  =  y2  (f)  +  q2  (f)  +  q2  (f) 

y:x  xy  xy^  xy^ 


(4. 3B) 


In  References  14  and  IS  the  multiple  coherence  function  for  a 
t hree- input/s ing! e  output  linear  model  is: 


y  (  f  )  - 

y:  x 


'  '  Vf)l  f1  •  *2yl(f)l  I’  •  y3y2!(f)1 


The  equivalence  of  these  latter  two  equations  is  straightforward 
to  verify. 

2 

I  lie  linear  coherence  function  y  (f)  in  Figure  IS  is  the  ordinary 

xy 

coherence  function  between  the  input  Z  ^ ( f )  =  X^ff)  and  the  total  output 
Y( f  )  ,  denoted  by : 


2  2 

y  (  n  =  y  ( n 

xy  z  y 


(4.40) 


The  nonlinear  coherence,  function  q  (f)  in  Figure  IS  can  be 

Xyb 

interpreted  as  the  ordinary  coherence  function  between  the  input  f) 
(  f)  and  the  total  output  Y(f),  denoted  by: 


q2:..  (f)  =  ..(f) 


V 


(4.41  ) 
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Similarly,  the  nonlinear  coherence  function  q  (f)  in  Figure  15 

yc 

can  he  interpreted  as  the  ordinary  coherence  function  between  the  input 
7,^(f)  =  |(f)  anr^  the  total  output  Yff),  denoted  by: 


2  2 
q  ...  (  f)  =  !T  ..(f) 


XY, 


V 


(4.42) 


Thus,  the  multiple  coherence  function  for  Equation  4.38  can  be 
written  as: 


y2  ff)  =  y2  (f)  +  y2  (f)  +  y2  (f)  (4.43) 

y :  X  7.  ^  y  Z  y  z  ^  y 

This  multiple  coherence  function  includes  all  terms  and  reflects 
the  goodness -of- f i t  of  the  postulated  model  at  each  frequency.  A 
cumulative  coherence  function,  defined  as  various  partial  sums  of  these 
component  coherences  for  a  select  number  of  terms,  is  presented  in  the 
next  section.  This  cumulative  coherence  allows  for  inspection  of  the 
model  for  select  combinations  of  terms. 

flood  models  will  be  obtained  at  desired  frequencies  when  the 
multiple  or  cumulative  coherence  function  is  close  to  unity.  The 
uncorrelated  output  noise  spectrum  G  (f)  representing  all  possible 
deviations  from  the  postulated  model  is  given  by  the  formula: 


G  ( f )  =  [1  -  y  ( f ) 1  G  (f) 

n n  y : x  yy 


(4.44) 


The  nonlinear  coherence  function  q  should  not  be  interpreted  as 

Xyb 

the  ordinary  coherence  function  between  the  input  X(f)  and  the  output 

Y^(f).  In  fact,  for  Gaussian  data  x(t),  the  ordinary  coherence  function 

between  X ( f )  and  Y^(f)  will  be  zero.  Similarly,  the  nonlinear  coherence 

function  q2  (  f)  should  not  be.  interpreted  as  the  ordinary  coherence 
xyc 

function  between  the  input  X ( f )  and  the  output  Y  (f).  Again,  for 
Gaussian  input  data,  this  ord'nary  coherence  function  is  also  zero. 
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4.5  Review  of  Reverse  Dynamic  System  MI/SO  Analysis  Technique 

Summarizing  to  this  point,  the  original  nonlinear  system  identifi¬ 
cation  problem  has  been  handled  as  follows: 

(a)  Perform  preliminary  analyses  to  detect  the  most  probable  non- 
linearities,  where  the  input  and  output  data  are  reversed. 

(b)  Postulate  a  nonlinear  differential  equation  of  motion.  It  is 
wise  at  this  stage  to  include  higher  order  terms  than  expected,  for  the 
following  reasons: 

1)  This  technique  will  correctly  determine  that  there  is  no 
frequency  response  function  for  terms  not  present  in  the  physical  model; 
calculating  this  zero-valued  function  increases  confidence  that  the  chosen 
order  of  the  model  accurately  represents  the  physical  system. 


2)  Inclusion  of  additional  terms  (e.g.,  adding  to  a  Duffing 
analysis)  does  not  bias  the  other  results;  in  other  words,  this  technique 
is  not  mode  1 -dependent . 

3)  The  only  cost  for  this  extra  information  is  slightly  in¬ 
creased  computer  analysis  costs  proportional  to  the  number  of  extra  input 
paths . 

(c)  Precompute  a  time  aeries  for  each  zero-memory  nonlinear 
function  gj(u). 

(d)  Compute  by  the  MI/SO  procedures  the  various  required  spectra 
for  solving  the  uncorrelated  form  of  the  model  as  shown  in  Figure  15. 


(e)  Compute  conditioned  frequency  response  functions  and  un¬ 
conditioned  frequency  response  functions  .  Pertinent  equations  are 
summarized  below;  the  frequency  dependence  of  all  functions  has  been 
deleted  for  clarity. 
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(4.45) 


I.,  =  G  /G 

1  V  z1*1 


T.„  =  G  /G 

2  z2y  V'2 


(4.46) 


T.,  =  G  /G 
3  Z3V  Z3Z3 


(4.47) 


Equations  4.26  through  4.28  relate  Lhese  h.  functions  to  the  A. 
1  S  J  J 


funct.  ions . 


( f )  Compute  component  and  total  spectral  functions  given  hy : 


i  ,  2 

g  =  r,,  f  g 

y  y  i  z  7. 

1  1 


(4.48) 


G  =  Lj  G 
y^v,  2' 


(4.49) 


G  =  T,  G 

y„y„  3  z.z. 


(4.50) 


G  =  G  +  G  +  G  +  G 
yy  yaya  ybyb  ycyc  nn 


(4.51) 


(g)  Compute  linear  and  nonlinear  coherence  functions  given  by: 


G  G 
ZjZ.j  yy 


(4.52) 


xy,  G 

Jb  yy 


G  G 
Z2Z2  yy 


(4.53) 


G  G 
Z3Z3  yy 


(4.54) 


The  last  three  equations  give  physical  meanings  to  the  three 
coherence  functions.  Equation  4.38,  the  sum  of  these  coherences,  is 
then  the  multiple  coherence  function  that  measures  the  goodness -of - f i t 
of  the  entire  model. 
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2  2 

The  reason  for  calling  q  and  q  nonlinear  coherence  functions 

xyb  Xyc  - 

is  because  the  output  spectrum  terms  G  and  G  are  the  result  of 

ybyb  ycyc 

the  original  input  X(f)  going  through  nonlinear  transformations  (e.g., 
squaring  and  cubing)  to  provide  the  correlated  inputs  and  X^  in 
Figures  13  and  14.  These,  in  turn,  lead  to  the  uncorrelated  inputs  Z^ 
and  Z^  in  Figure  15  from  which  one  calculates  the  uncorrelated  outputs 
Y,  and  Y  to  obtain  G  and  G 

b  c  Vb  Vc 

2 

The  first  cumulative  coherence  function  is  defined  as  IT  (f).  This 

xy 

identifies  the  proportion  of  the  output  spectrum  due  to  x(t)  passing 

through  the  optimum  linear  system  L^(f).  The  second  cumulative 

2  2 

coherence  function  is  defined  by  Y  (f)  +  q  (f).  This  identifies  the 

- - -  —  xyv  xyb 

proportion  of  the  output  spectrum  due  to  x(t)  and  z^Ct)  passing  through 

L^(f)  and  L^Cf).  The  third  cumulative  coherence  function  defined  by 
2  2  2 

Y  (f)  +  q  (f)  +  q  (f)  is  the  same  here  as  the  multiple  coherence 


function. 


25 


5.0  NUMERICAL  RESULTS  OF  REVERSE  DYNAMIC  SYSTEM  MI/SO  ANALYSIS 

5 . 1  Overview 

An  evaluation  of  the  MI/SO  method  applied  to  reverse  nonlinear 
dynamic  systems  has  been  conducted  for  the  SDOF  systems  described  in 
Section  2.  The  evaluation  is  presented  here  in  three  sections,  namely: 

(a)  All  systems  excited  by  a  random  force  with  unit  variance, 
without  corruption  due  to  noise,  to  demonstrate  applicability  of  the 
method . 

(b)  The  Duffiiig-Van  der  Pol  system  excited  by  various  levels  of 
RMS  excitation,  without  corruption  due  to  noise,  to  verify  if  the  results 
are  Independent  of  excitation  energy. 

(c)  The  Duffing-Van  der  Pol  system  excited  by  a  random  force  with 
unit  variance  corrupted  by  various  levels  of  noise,  to  evaluate  the 
effects  of  extraneous  noise. 


The  parameters  used  in  the  numerical  simulations  are  listed  in 
Table  1. 

Table  1.  Simulated  System  Parameters* 


System  Type 

k 

c 

f 

n 

C 

Nonlinear 

Parameters 

Linear 

355.3 

3.77 

3.0 

0.  10 

- 

- 

Duff ing 

355.3 

3.77 

3.0 

0.  10 

k3=2  x  107 

c3=0 

Van  der  Pol 

355.3 

3.77 

3.0 

0.  10 

k3=o 

c3=2  x 

io6 

Combined  D-VdP 

355.3 

3.77 

3.0 

0.  10 

k3=2  x  107 

c3=2  x 

io6 

Mathieu 

355.3 

3.77 

3.0 

0.  10 

k  =177.7 
m 

f  =1.0 
m 

Dead-Band 

355.3 

1.89 

3.0 

0.05 

u  =0.0025 

g 

Extreme  D-VdP 

0 

0 

0 

0 

k3=2  x  107 

c3=2  x 

io6 

*A11  systems  used  unit  mass  m=1.0. 
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5  2  Identification  of  Nonlinear  Systems 


This  first  series  of  simulations  was  selected  to  demonstrate  that 
the  MI/SO  method  is  capable  of  correctly  identifying  the  linear  and  non¬ 
linear  properties  of  a  wide  variety  of  nonlinear  dynamic  systems  (with 
ccnctauL  coetr icients ) .  The  numerical  techniques  used  in  these  simula¬ 
tions  were  described  in  Sections  2.1  and  2.2.  The  applied  excitation 
force  used  in  all  cases  is  illustrated  in  Figures  16  and  17,  which  in¬ 
cludes  a  time  history  segment,  the  force  autospectrum,  and  the  force 
histogram.  The  total  time  history  excitation  and  response  for  each 
simulation  is  8,192  samples  (163.84  seconds),  with  subrecord  lengths  of 
512  for  the  Fourier  transforms  (10.24  seconds).  Therefore,  each  spec¬ 
tral  and  frequency  response  function  is  based  on  16  averages  with  a 
bandwidth  resolution  of  0.0977  Hz.  As  is  common  practice  in  linear 
spectral  analysis,  a  Hanning  window  was  used  to  reduce  the  effects  of 
leakage;  a  256-point  overlap  was  employed  to  recover  information  "lost" 
due  to  the  windowing. 

Table  2  cross- lists  the  graphical  results  presented  in  the  sub¬ 
sequent  sections  according  to  the  type  of  function  and  type  of  nonlinear 
system.  This  table  is  presented  to  allow  for  quick  inspection  of  the 
results . 


Table  2.  Comparison  of  Linear  and  Nonlinear  Results 


Figure  Numbers 

System  Type 

Reverse  Frequency 
Response  Functions: 

Cumulative 

Coherences 

Response 

Histogram 

Linear 

Nonlinear 

SDOF  Linear 

18a, 

18b 

- 

19 

20 

Duff ing 

21a, 

21b 

22a,  22b 

23 

24 

Van  der  Pol 

25a, 

25b 

26a,  26b 

27 

28 

Duffing-Van  der  Pol 

29a, 

29b 

30a,  30b 

31 

32 

Mathieu 

33a, 

33b 

34a,  34b 

35 

36 

Dead-Band 

Extreme  Duffing- 

37a, 

37b 

38a,  38b 

39 

40 

Van  der  Pol 

41a, 

41b 

42a,  42b 

43 
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5.2.1  SDOF  Linear  System.  Results  are  presented  first  for  a  SDOF 
linear  system  to  illustrate  the  new  technique  for  this  familiar  system 
type.  The  linear  reverse  dynamic  system  frequency  responses,  magni¬ 
tude  and  phase,  are  presented  in  Figure  18a  and  18b  where  f  =3.0  Hz  and 

n 

C  =  0.10.  This  function  has  been  calculated  by  the  usual  linear  SI/SO 
technique  (dashed  line)  and  the  new  technique  (solid  line)  which  assumed 
a  simple  Duffing  model.  As  one  would  expect,  there  is  virtually  no  dif¬ 
ference  between  results  from  the  two  methods  since  there  is  no  nonlinearity 
in  the  system.  The  cumulative  coherence  function  provided  in  Figure  19, 
which  is  near  unity,  confirms  the  fact  that  the  system  is  linear  with  no 
significant  noise  corruption.  The  dip  in  coherence  near  resonance  is 
due  to  the  well-known  effect  of  leakage.  The  response  histogram  provided 
in  Figure  20  exhibits  an  apparent  Gaussian  distribution  which  further 
confirms  the  conclusion  that  the  system  was  linear. 

The  linear  frequency  response  function  A^  is  the  reciprocal  of  the 
usual  SI/SO  function  H.  The  SDOF  system  stiffness,  k  (exact  value  = 

355.5),  is  obtained  from  the  asymptotic  behavior  of  A^  for  f  — >  0.  In 
addition,  the  natural  frequency,  f^  (exact  value  =3),  is  estimated  based 

on  the  local  minimum  value  of  A  .  Consequently,  the  system  mass,  m,  is 

2  1 

estimated  using  k/ ( 2tt f ^)  .  The  value  of  the  viscous  damping  is  estimated 

on  the  basis  of  the  magnitude  of  A,  at  f=f  where  I  A,  I  =  2-irf  c. 

1  n  1  1  1  u 

5.2.2  Duffing  SDOF  Nonlinear  System.  Results  for  a  Duffing  SDOF 

system  with  identical  linear  system  parameters  as  the  above  linear  system, 

and  cubic  stiffness  k^  =  2  x  10^,  are  presented  next.  For  this  analysis, 

two  inputs  are  required  (as  discussed  in  Section  4.1):  u(t),  and  the 

3 

cubic  displacement  u  (t).  Estimates  of  the  linear  frequency  response  A^ 
magnitude  and  phase  are  provided  in  Figure  21a  and  21b.  The  results 
(solid  curve)  correctly  identify  this  function  (since  they  agree  with 
the  results  in  Figure  18a  and  18b  from  the  previous  section).  The  SI/SO 
based  frequency  response  function  returns  a  "best  linear  fit"  that  is 
corrupted  by  the  correlated  cubic  stiffness  parameter;  inspection  of 
that  function  shows  that  it  incorrectly  estimates  k  =  500  (exact=355) 
and  f  =  3.6  (exact=3).  Because  it  assumes  a  linear  system  model,  the 
usual  SI/SO  results  are  not  capable  of  correctly  identifying  parameters 


28 


that  are  invariant  with  excitation  variance;  in  fact,  estimates  of  k  and 
m  are  incorrect  for  all  input  variances.  While  this  particular  SI/SO 
frequency  response  function  would  hold  for  all  realizations  with  the 
same  variance,  it  would  vary  for  different  excitation  levels,  and  is 
therefore  an  invalid  model  for  this  system. 

The  MI/SO  method,  however,  identifies  the  nonlinear  component  of 
the  reverse  system  frequency  response,  A^,  as  shown  in  Figure  22a  and 
22b  in  terms  of  its  real  and  imaginary  parts.  This  function  is  correctly 
identified;  the  real  part  is  equal  to  2  x  10^  (as  used  in  the  simulations), 
and  the  imaginary  part  is  essentially  zero.  Thus,  the  MI/SO  technique 
correctly  identified  this  system  as  a  Duffing  system.  The  validity  of 
these  results  is  confirmed  through  inspection  of  the  cumulative  coherence 
function  shown  in  Figure  23,  which  is  nearly  unity  over  the  entire  fre¬ 
quency  range.  Furthermore,  Figure  23  shows  that  approximately  90  percent 
of  the  output  passes  through  the  strictly  linear  path  given  by  L^.  The 
fact  that  this  system  description  is  excitation-invariant  is  examined  in 
Section  5.3. 

The  displacement  histogram  for  the  Duffing  response  relative  to  the 
unit  input  variance  is  provided  in  Figure  24.  The  distortion  from  Gaussian 
is  rather  slight,  so  it  would  not  serve  well  as  a  preliminary  analysis 
tool  for  identification  of  this  type  of  nonlinearity. 

5.2.3  Van  der  Pol  SDOF  Nonlinear  System.  The  basic  behavior  of  a 

Van  der  Pol  SDOF  system  with  linear  system  parameters  identical  with  the 

system  described  in  Section  5.2.1  and  cubic  damping  with  c  =  2  x  10^  is 

3  3 

presently  examined.  As  with  the  Duffing,  two  inputs  (u  and  u  )  are  used. 
Estimates  of  the  linear  reverse  dynamic  system  frequency  response  function 
magnitude  and  phase  are  provided  in  Figure  25a  and  25b.  The  MI/SO  estimate 
of  A^  (solid  curve)  correctly  identifies  this  function,  while  the  SI/SO 
estimate  (dashed  curve)  again  identifies  an  effective  linear  system  cor¬ 
rupted  by  the  nonlinear  damping  component.  As  with  the  SI/SO  result  for 
the  Duffing  system,  this  estimate  is  not  invariant  with  input  variance. 

The  MI/SO  technique  correctly  identifies  the  nonlinear  component  A^ 
as  shown  in  Figure  26a  and  26b  in  terms  of  the  real  and  imaginary  com- 

3 

ponents.  Given  the  nonlinear  input  of  u  ,  and  the  fact  that  the  Van  der 
Pol  nonlinearity  is  proportional  to  the  out-of-phase  velocity,  the 
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expected  frequency  response  function  would  have  a  negligable  real  component 
and  an  imaginary  component  equal  to  2  x  10^  times  ( 2tt f / 3 ) .  Inspection  of 
the  log  scale  in  Figure  26b  shows  that  the  MI/SO  estimates  are  correct. 

Quantitative  confirmation  of  the  estimates  is  confirmed  via  the 
cumulative  coherence  function  shown  in  Figure  27.  The  lower-valued  curve 
indicates  that  nearly  all  of  the  excitation  for  this  system  below  2  Hz 
is  associated  with  the  linear  path  L^.  Above  2  Hz,  the  contribution 
from  the  nonlinear  term  increases.  The  validity  of  the  whole  model  is 
confirmed  since  the  cumulative  coherence  is  close  to  unity  over  all 
frequencies.  The  displacement  histogram  for  this  response  is  shovfn  in 
Figure  28,  showing  strong  deviations  from  the  Gaussian  shape  expected 
for  a  linear  system. 

5.2.4  Combined  Ihif f ing-Van  der  Pol  SDOF  Nonlinear  System.  The 

linear  and  nonlinear  parameters  used  here  were  identical  to  the  previous 
simulations.  Estimates  of  the  linear  reverse  frequency  response  func¬ 
tions  are  shown  in  Figure  29a  and  29b.  As  expected,  the  SI/SO  estimate 
is  contaminated  by  the  nonlinear  paths,  while  the  MI/SO  estimate  is  ac¬ 
curate  . 

The  estimate  for  the  nonlinear  reverse  frequency  response  function 
A^  is  provided  in  Figure  30a  and  30b  in  terms  of  the  real  and  imaginary 
functions.  Comparison  to  the  previous  Duffing  and  Van  der  Pol  values 
shows  that  the  parameters  are  again  properly  identified.  Inspection  of 
the  cumulative  coherence  shown  in  Figure  31  confirms  this  conclusion. 

Note  that  the  histogram  for  the  displacement  response  shown  in 
Figure  32  appears  close  to  Gaussian,  even  though  there  are  two  addi¬ 
tional  nonlinear  paths  for  this  system.  Therefore,  this  would  not 
provide  any  useful  information  for  evaluating  whether  the  system  was 
linear  or  not;  indeed,  it  would  probably  reinforce  the  erroneous  con¬ 
clusion  that  the  system  was  linear. 

An  important  point  can  be  made  here  on  the  value  of  the  MI/SO  method 
for  nonlinear  system  identification  problems.  Even  though  the  absolute 
contribution  from  the.  nonlinear  terms  is  small  for  this  particular  simu¬ 
lation  (as  shown  by  the  coherence  functions)  and  the  response  is  therefore 
nearly  linear,  the  MI/SO  technique  successfully  identified  the  system  as 
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nonlinear.  Thus,  the  method  appears  to  be  a  robust  tool  of  great  value. 
Conversely,  a  standard  SI/SO  analysis  of  the  coherence  function  would 
almost  assuredly  have  concluded  that  the  system  behaved  in  a  linear 
manner;  this  would  clearly  lead  to  erroneous  results  for  higher  values 
of  excitation  variance. 

5.2.5  Mathieu  SDOF  Nonlinear  System.  Nonlinear  systems  described 
by  the  Mathieu  equation  exhibit  temporal  rather  than  amplitude-dependent 
nonlinear  behavior;  the  fishtailing  behavior  of  vessels  at  a  single  point 
mooring  is  a  good  example. 

The  presence  of  periodic  variation  in  stiffness  is  first  detected 
by  preliminary  data  analysis  as  discussed  in  Sections  3.1  and  3.2.  This 
information  is  needed  to  apply  the  new  technique. 

Employing  the  estimate  for  f  ,  three  reverse  dynamic  inputs  are 
computed  for  the  Mathieu  analysis  as  specified  in  Equation  4.9.  The 
three  frequency  response  functions  A^,  A^,,  and  A^  are  then  identified 
using  the  new  method.  For  these  simulations,  the  choice  of  parameters 
resulted  in  negligable  output  from  the  A3  system.  In  an  actual  Mathieu 
system,  the  phase  of  the  time-varying  stiffness  is  unknown  (it  is 
directly  related  to  the  arbitrary  start  time  for  collecting  data);  hence, 
both  a  cosine  and  sine  component  are  modeled  according  to  Equation  4.9. 
For  these  simulations,  =  0  was  used,  removing  the  odd  component  pro¬ 
portional  to  the  sine,  and  thereby  effectively  eliminating  the  A3  path. 

Estimates  of  the  reverse  dynamic  linear  frequency  response  function, 
Aj  magnitude  and  phase,  are  provided  in  Figure  33a  and  33b.  It  is  note¬ 
worthy  that  the  estimate  based  on  the  SI/SO  technique  (dashed  curve)  is 
quite  erroneous,  indicating  apparent  additional  resonant  frequencies  at 
the  sum  and  difference  frequencies.  However,  the  coherence  associated 
with  the  SI/SO  estimate  (not  shown  here)  is  quite  low,  indicating  that 
the  model  is  probably  not  linear.  Once  again,  the  MI/SO  linear  frequency 
response  function  is  correct. 

The  MI/SO  estimate  for  A^  is  shown  in  Figure  34a  and  34b  in  terms 
of  the  real  and  imaginary  parts.  The  correct  values  for  these  functions 
are  seen  to  be  RetA^]  =  178  and  ImlA^]  =  0  from  Table  1  and  Equation 
4.9;  the  estimate  is  excellent  over  the  portion  of  the  frequency  range 
with  the  predominant  energy. 
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As  before,  quantitative  confirmation  is  also  provided  by  the  cumu¬ 
lative  coherence  functions  shown  in  Figure  35.  First,  the  cumulative 
coherence  is  essentially  unity  over  the  whole  range,  validating  the 
Mathieu  model.  Second,  the  linear  component  of  the  response  that  passes 
through  Lj  is  rather  low,  indicating  that  a  strictly  linear  model  would 
not  provide  a  very  useful  model  for  a  Mathieu  system. 

The  displacement  histogram  for  the  Mathieu  response  is  provided  in 
Figure  36.  The  somewhat  flattened  top  section  is  indicative  that  the 
response  may  have  some  nonlinear  contributions. 

5.2.6  Dead-Band  SDOF  Nonlinear  System.  Some  dynamic  systems  (for 
example,  a  loosely-moored  ship)  only  exhibit  a  stiffness  when  the  res¬ 
ponse  reaches  a  certain  value.  The  identification  and  characterization 
of  such  systems  is  explored  here. 

A  SDOF  system  with  a  dead-band  spring  may  be  crudely  modeled  as  a 
Duffing  system.  If  MI/SO  analysis  is  performed  on  dead-band  responses 
with  an  assumed  Duffing  model,  a  cubic  stiffness  will  be  identified. 

Even  though  the  cumulative  coherence  function  associated  with  this 
assumed  system  would  show  dramatic  improvement  over  the  best  linear  fit, 
it  will  never  approach  unity  if  the  dead  band  is  of  the  same  order  of 
magnitude  as  the  RMS  responses.  Moreover,  the  identified  Duffing  system 
will  not  be  Invariant  to  the  input  excitation.  It  is  imperative  to  cor¬ 
rectly  identify  the  system  as  having  a  dead  band  and  to  incorporate  that 
information  into  the  MI/SO  model. 

Assume  that  preliminary  data  analysis  has  correctly  identified  the 
value  of  u^  that  defines  the  dead  band.  A  MI/SO  model  based  on  Equation 
4.10  can  then  be  formulated  (with  two  inputs). 

Estimates  of  the  reverse  linear  frequency  response  function,  A^ 
magnitude  and  phase,  are  provided  in  Figure  37a  and  37b,  where  the  dashed 
curve  shows  the  best-fit  linear  SI/SO  estimate  and  the  MI/SO  estimate  is 
shown  as  the  solid  curve.  The  SI/SO  estimated  natural  frequency  is  too 
low,  which  is  expected  since  the  best  linear  stiffness  averaged  over  the 
entire  response  range  will  be  very  low.  The  MI/SO  estimate  is  correct. 

Figure  38a  and  38b  show  the  real  and  imaginary  components  of  A^.  A 
relatively  accurate  estimate  of  the  linear  stiffness  value  (-355)  is 
obtained  for  all  frequencies,  with  a  maximum  error  of  7  percent  in  the 
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2-  to  3-Hz  band.  The  cumulative  coherence  functions  for  the  MI/SO 
estimates  are  shown  in  Figure  39.  The  linear  value  of  approximately  65 
percent  approaches  unity  when  the  nonlinear  contribution  is  added  (the 
somewhat  lower  coherence  near  resonance  is  due  to  leakage  because  of  low 
damping  in  the  system) .  The  response  histogram  is  provided  in  Figure 
40,  showing  markedly  non-Gaussian  form. 

5.2.7  Extreme  Duffing-Van  der  Pol  SDOF  Nonlinear  System.  All  of 

the  nonlinear  systems  studied  up  to  this  point  were  constructed  by  ap¬ 
pending  nonlinear  terms  onto  the  original,  linear,  differential  equation 
of  motion.  And,  as  shown  by  the  various  coherence  functions  presented 
for  each  case,  the  linear  contribution  typically  dominated  the  total 
nonlinear  response.  This  is  not  unexpected  for  the  types  of  simple  non¬ 
linear  systems  examined.  This  final  simulation  was  constructed  as  an 
extremely  nonlinear  system  to  more  fully  demonstrate  the  capabilities  of 
the  MI/SO  technique  compared  to  usual  SI/SO  procedures.  As  stated  in 
Section  2.1,  this  final  system  consists  of  the  inertial  term,  with  cubic 
stiffness  and  damping  only  (as  given  by  the  combined  Duffing-Van  der  Pol 
system  with  the  linear  stiffness  and  damping  terms  zeroed  out). 

Since  the  system  lacks  a  linear  stiffness,  the  natural  frequency 
(see  Equation  4.4)  is  identically  zero.  However,  for  a  given  level  of 
excitation  variance,  this  cubic  stiffness  would  appear  as  a  linear  stiff¬ 
ness  to  a  best-fit  linear  analysis;  thus,  a  usual  SI/SO  analysis  would 
erroneously  estimate  an  apparent  (but  false)  natural  frequency.  In  this 
case,  the  SI/SO  results  would  lead  to  serious  quantitative  errors  (the 
system  mass  and  linear  stiffness)  as  well  as  qualitative  errors  (the 
basic  system  behavior). 

Estimates  of  the  reverse  linear  frequency  response  functions,  mag¬ 
nitude  and  phase,  are  shown  in  Figure  41a  and  41b.  The  MI/SO  estimate 

2 

(solid  curve)  correctly  identifies  the  linear  system  as  -(2nf)  m.  The 
SI/SO  estimate  (dashed  curve)  is  completely  wrong  for  the  reasons  dis¬ 
missed  above.  The  MI/SO  estimate  for  A^,  shown  in  Figure  42a  and  42b, 
is  correct  as  seen  by  comparison  to  previous  cases.  The  goodness-of-f it 
of  the  MI/SO  model  is  confirmed  by  the  cumulative  coherence  function 
shown  in  Figure  43. 
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5.3  Results  from  Different  Excitation  Levels 


It  was  stated  in  Section  5.1  that  the  reverse  frequency  response 
functions  identified  using  the  MI/SO  method  are  globally  valid,  that  is, 
that  the  1 s  are  independent  of  the  input  variance.  This  is  a  required 
characteristic  for  any  valid  nonlinear  model.  Moreover,  it  was  stated 
that  the  effective  linear  system  estimates  deduced  from  SI/SO  analysis 
are  valid  only  for  the  particular  excitation  level  encountered.  These 
statements  are  verified  in  this  section  by  simulating  the  combined 
Duffing-Van  der  Pol  nonlinear  case  (see  Section  5.2.4)  with  a  range  of 
excitation  variances. 

F jr  these  simulations,  the  combined  system  is  subjected  to  an 
excitation  with  mean  square  values  of  0.2,  0.5,  1.0,  and  2.0  (the  unity 
variance  results  were  previously  presented  in  Section  5.2.4).  The 
frequency  response  functions  and  k^  and  the  cumulative  coherence 
functions  are  provided  for  these  cases  as  follows: 


Input  variance: 

0.2 

0.5 

1.0 

2. 

.0 

Figure  Numbers: 

44-46 

47-49 

29-31 

50- 

-52 

The  estimates  of  A^,  in  order  of  increasing  force  variance,  are 
shown  in  Figures  44,  47,  29,  and  50,  respectively,  for  the  MI/SO  (solid 
curve)  and  SI/SO  (dashed  curve)  methods.  The  MI/SO  estimate  is  invariant 
while  the  SI/SO  estimates  erroneously  exhibit  an  increasing  apparent 
resonance  (as  expected).  Note  that  both  estimates  are  similar  for  the 
lowest  force  level  where  the  nonlinear  effects  are  minimal. 

The  MI/SO  estimates  of  A^,  again  in  order  of  increasing  force 
variance,  are  shown  in  Figures  45,  48,  30,  and  51.  These  results  are 
also  invariant  with  force  level.  The  cumulative  coherences  for  all  these 
simulations  are  shown  in  Figures  46,  49,  31,  and  52.  The  near-unity 
values  confirm  the  validity  of  the  results,  and  in  turn,  the  validity  of 
the  MI/SO  identification  method  for  these  cases. 
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5.4  Results  from  Different  Extraneous  Noise  Levels 


A  practical  system  identification  method  must  be  capable  of  detect¬ 
ing  signals  in  the  presence  of  noise.  The  purpose  of  the  simulations  in 
this  section  is  to  numerically  evaluate  the  robustness  of  the  MI/SO  method 
to  various  levels  of  noise. 

Extraneous  noise  may  occur  in  both  input  and  output  measurements. 

In  the  case  of  linear  systems  subjected  to  multiple  excitations,  the 
noise  source  is  usually  assumed  at  the  output  (response)  signal  only. 

Such  a  noise  type  produces  random  errors  and  no  biasing  errors  (Ref  15). 

In  order  to  be  consistent  with  this  assumption,  the  present  study  assumes 
a  noise  source  in  the  reverse  dynamic  system  inputs  u^(t),  which  correspond 
physically  to  the  displacement  response.  This  type  of  noise,  which  does 
not  pass  through  the  system,  theoretically  produces  random  errors  only. 

For  the  purpose  of  investigating  the  effects  of  noise  contamination, 
the  combined  Duffing-Van  der  Pol  system  (see  Section  5.2.4)  was  chosen. 
Noise-free  estimates  of  A.  and  the  cumulative  coherence  were  presented 
in  Figures  29  through  31.  Analyses  of  this  system  with  mean  square  noise 
levels  versus  displacement  response  of  5,  10,  and  20  percent  are  presented. 
The  numerical  noise  was  approximately  "white"  as  discussed  in  Section 
2.  1. 

The  frequency  response  function  estimates  A^  and  A^  and  cumulative 
coherences  for  these  cases  are  provided  as  follows: 


Noise  Level  (%) 
Figure  Numbers 


0 

5 

10 

20 

29-31 

53-55 

56-58 

59-61 

Examination  of  the  A^  estimates  indicates  that,  although  random 
errors  increase  with  noise  level,  biasing  errors  are  not  significant. 

The  A^  estimates  indicate  an  increasing  level  of  both  random  and  biasing 
errors,  especially  in  the.  imaginary  component,  as  the  noise  level  in¬ 
creases.  The  cumulative  coherence  functions,  while  decreasing  slightly 
as  the  noise  increases,  still  provide  a  quantitative  confirmation  that 
the  model  is  accurately  identifying  the  system  properties. 
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6 . 0  CONCLUS IONS 


This  report  has  addressed  the  problem  of  identifying  the  properties 
of  nonlinear  systems  from  measurements  of  input  and  output  stochastic 
data.  This  new  method  for  nonlinear  system  identification  is  based  upon 
mul t ipl e- input/s ingl e-output  linear  procedures  applied  to  reverse  dynamic 
systems  that  are  defined  for  proposed  nonlinear  differential  equations 
of  motion.  It  is  shown  that  this  new  method  allows  for  efficient  prac¬ 
tical  identification  of  a  wide  variety  of  nonlinear  system  behavior. 

Conclusions  from  this  investigation  regarding  the  characteristics 
and  .  of  this  method  are: 

1.  The  method  converts  the  original,  single-input/single-output 
nonlinear  problem  into  an  equivalent  multiple-input/single-output  linear 
problem,  thereby  allowing  for  the  use  of  standard  linear  theory  and  es¬ 
tablished  computer  programs.  Nonlinear  system  amplitude-domain  proper¬ 
ties  are  determined  separate  from  nonlinear  system  frequency-domain 
properties . 

2.  The  method  computes  linear  frequency  response  functions  for  each 
nonlinear  term  in  the  proposed  nonlinear  differential  equation  of  motion. 
The  relative  contributions  from  each  individual  term  as  well  as  the  en¬ 
tire  postulated  model  can  be  quantitatively  assessed  through  appropriate 
cumul''tive  coherence  functions. 

3.  The  method  identifies  specific  terms  as  postulated  in  the 
proposed  nonlinear  differential  equation  of  motion  However,  since  any 
number  of  terms  can  be  investigated,  enough  terms  can  be  included  to 
assure  convergence,  and  the  results  need  not  be  model -dependent . 

The  method  does  not  require  any  assumptions  regarding  the 
probability  distributions  or  spectral  shapes  of  the  input  or  output 
data . 
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5.  The  model  gives  results  that  are  independent  of  the  actual 
input  excitation  level  as  required  for  valid  nonlinear  models.  Also, 
the  results  are  not  sensitive  to  typical  noise  levels. 

6.  The  method  correctly  identifies  nonlinear  system  properties 
even  for  low  levels  of  excitation  where  conventional  single-input/ 
single-output  linear  analysis  would  show  a  high  ordinary  coherence  and 
erroneously  lead  to  the  conclusion  that  the  system  is  linear.  Thus,  the 
new  method  improves  system  identification  even  for  nearly-linear  systems. 

7.  Various  types  of  nonlinear  system  differential  equations  of 
motion  were  examined,  including  Duffing,  Van  der  Pol,  Mathieu,  and 
Dead-Band.  All  were  successfully  identified  by  applying  this  method  to 
their  reverse  dynamic  systems. 
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7.0  RECOMMENDATIONS 


1.  Investigate  force  excitations  with  arbitrary  spectral  shapes 
(such  as  the  Pierson-Moskowitz  wind-wave  model)  and  apply  these 
excitations  to  various  nonlinear  systems. 

2.  Simulate  extended  SDOF  nonlinear  systems  with  more  complicated 
nonlinear  restoring  and  damping  forces  than  those  considered  here. 

3.  Study  effects  of  noise  in  greater  detail  for  cases  of  noise  at 
both  the  inputs  and  outputs  of  simulated  SDOF  systems. 

4.  Develop  procedures  for  analyzing  and  identifying  different  terms 
in  MDOF  systems  and  associated  MDOF  nonlinear  differential  equations  of 
motion . 

5.  Demonstrate  these  procedures  by  investigating  the  properties  of 
an  actual  nonlinear  physical  system  using  measured  data. 

6.  Investigate  the  use  of  this  new  method  for  the  nonlinear  system 
identification  of  nonlinear  equations  of  motion  where  the  coefficients 
are  not  necessarily  constants. 
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FORCE 


SDOF  Nonlinear  Equations 


Mu(t)  +  Cu( t)  +  Ku( t)  +  P(u,u,t)  -  F(t) 

F  r(  u ,  u ,  t ) 

where 

Fr(u,u,t)  -  F(t)  -  MU(t) 

Figure  1.  SDOF  nonlinear  system  and  equations. 


LINEAR  FORCE  VS.  0*SPt_ 


RESTORING  FORCE  VS  DISPLACEMENT  LINEAR  SYSTEN 


WSPL 


*10-* 


Figure  2.  Restoring  force/displacement  patterns  for  Linear  system. 
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DUFFING  FORCE  VS.  DISPL. 


RESTORING  FORCE  VS.  DISPLACEMENT.  .CUFFING  SYSTEM 


Figure  5.  Restoring  force/displ aoement  patterns  for  combined 
Duffing-Van  der  Pol  system. 
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mathieu  force  vs.  oispl 


RESTORING  FORCE  VS.  DISPLACEMENT.. MATHIEU  SYSTEM 


Figure  6.  Restoring  force/displacement  patterns  for  Mathieu  system. 


CAPPING  FORCE  VS.  DISPL. 


RESTORING  FORCE  VS.  DISPLACEMENT..  GAP  PING  SYSTEM 


S 


-0.015  -0.01  -O.OOS  O  0.005  0.01  0.015 
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Figure  7.  Restoring  force/displacement  patterns  for  Dead-Band  system. 


AUTOSPECTRUM  OP  RESPONSE.. .MATHIEU  SYSTEM 


§  io-4 


0  0.5  1  I.S  2  2-5  3  3.5  4  4.5  3 

TWO  (HI) 

Figure  8.  Autospectrum  for  Mathieu 
system  response. 


p(u,ii,t) 


■u(t)  +•  cu(t)  ♦  ku(t)  +■  p(u,u,c)  -  F(t) 


vh*c*:  F(e)  -  fore*  Input 

u(c)  -  dliplccaaant  output 

Figure  9.  Single-input/single-output 

nonlinear  model  with  feedback 
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u(C) 


Nonlinear 


System 


p(u,u,t) 


mil(t)  +  cu(t)  +  ku(t)  +  p(u,u,t)  -  F(t) 


gj^uJ-uCt) 


V)— ►F(t) 


where:  u(c)  -  displacement  input  gj(u) - ►  Aj(f)  - ' 

F(t)  -  force  output 

Figure  10.  Reverse  single-input/single-  Figure  11.  Three- input/s ingle-output 

output  nonlinear  model  without  linear  model  where  g2(u)  and 

feedback.  g^(u)  are  zero*memory  nonlinear 

functions  of  u. 


Squarer  - ke 

-  Xjif)  - 


Figure  12.  SI/SO  nonlinear  model  where 

nonlinear  stiffness  terms  have 
constant  coefficients. 


Figure  13.  Three-input/single-output 

linear  model  where  nonlinear 
stiffness  terms  have  constant 
coefficients . 


Z2(f>  -  XjU) 


z2(f>  -  Xj.i(f)- 
z3(f)  -  x3.2,(f) 


Figure  14.  Three-input/single-output 

linear  model  where  nonlinear 
stiffness  terms  have  frequency- 
dependent  coefficients. 


Figure  15.  Equivalent  three-input/single¬ 
output  linear  model  to  Figure 
14  with  uncorrelated  inputs. 


COH  ABS(HAY) 


Figure  18.  SDOF  linear  system  impedance 
(magnitude  and  phase). 
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Figure  19. 


SDOF  linear  system  coherence 
function . 


-drfl 


Ebb. 


-0  -0  -4  -2  0  2  4  8 

DISPLACEMENT 

Figure  20.  SDOF  linear  system 

displacement  histogram. 


COH 


LINEAR  IMPEDANCE  FUNCTION  MAGNITUDE  *  SIMPLE  LINEAR  FRF  EST. 


UNEAR  IMPEDANCE  FUNCTION  PHASE  k  SIMPLE  UNEAR  FRF  EST. 


Figure  21.  Duffing  system  linear  impedance 
(magnitude  and  phase). 


Il07  REAL  PART  OF  NONUNEAR  IMPEDANCE  FUNCTION 


1 10*  IMAGINARY  PART  OF  NONUNEAR  IMPEDANCE  FUNCTION 


Figure  22.  Duffing  system  nonlinear 

impedance  (real  and  imaginary). 


CUMULATIVE  COHERENCE  FUNCTION  DISPL  HISTOGRAM  STD  =  0  001  8102  SAMPLES 


FREQ  (HZ) 

Figure  23.  Duffing  system  cumulative 
coherence  functions. 


DISPLACEMENT  Il0> 

Figure  24.  Duffing  system  displacement 
histogram. 


48 


COH  BEAHHBY)  ABS(HAY) 


IJNEAR  IMPEDANCE  FUNCTION  PHASE  *  SIMPLE  LINEAR  FRF  E'  T 


COH  REAUHBY)  ABS(HAY) 


LINEAR  IMPEDANCE  FUNCTION  MAGNITUDE  it  SIMPLE  LINEAR  FRF  EST 


UNEAR  IMPEDANCE  FUNCTION  PHASE  It  SIMPLE  UNEAR  PRF  EST 


Figure  29.  Combined  Duffing-Van  der  Pol  system  linear  impedance 
(magnitude  and  phase). 


Iio7  REAL  PART  of  nonunear  impedance  function 


1 10’  IMAGINARY  PART  OP  N0NUNEAR  IMPEDANCE  FUNCTION 


Figure  30.  Combined  Duffing-Van  der  Pol  system  nonlinear  impedance 
(real  and  imaginary). 


CUMULATIVE  COHERENCE  FUNCTION 


Figure  31.  Combined  Duffing-Van  der  Pol 
system  cumulative  coherence 
functions . 


DISPL  HISTOCRAM  STD  =  0  002  8192  SAMPLES 


DISPLACEMENT 

Figure  32.  Combined  Duffing-Van  der  Pol 
system  displacement  histogram. 
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REAUHBY)  ABS(HAY) 


UNEAR  IMPEDANCE  FUNCTION  MAGNITUDE  it  SIMPLE  LINEAR  FRF  EST 


UNEAR  IMPEDANCE  FUNCTION  PHASE  fc  SIMPLE  UNEAR  FRF  EST 
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Figure  33.  Mathieu  system  linear  impedance 
(magnitude  and  phase) . 
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Figure  34. 
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Mathieu  system  nonlinear 
impedance  (real  and  imaginary). 
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Figure  35. 
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Mathieu  system  cumulative 
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Figure  36.  Mathieu  system  displacement 
histogram. 
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Figure  45.  Nonlinear  impedance  real  and 
imaginary  for  0.2  input 
variance . 
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Figure  48.  Nonlinear  impedance  real  and  imaginary 
for  0.5  input  variance. 
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Figure  50.  Linear  impedance  magnitude  for 
2.0  input  variance. 
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Figure  51.  Nonlinear  impedance  real  and 
imaginary  for  2.0  input 
variance. 
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Figure  52. 
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Cumulative  coherence  functions 
for  2.0  input  variance. 


Figure  53. 


FREQ  (HZ) 

Linear  impedance  magnitude 
for  5%  noise. 
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Figure  54. 
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Nonlinear  impedance  real  and 
imaginary  for  5%  noise. 
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Figure  55.  Cumulative  coherence 

functions  for  5%  noise. 
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Figure  56.  Linear  impedance  magnitude 
for  10%  noise. 
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Figure  57. 

Nonlinear  impedance  real  and 
imaginary  for  10%  noise. 
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Figure  58.  Cumulative  coherence  functions 
for  10%  noise. 
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Figure  59.  Linear  impedance  magnitude 
for  20%  noise. 
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Figure  60.  Nonlinear  impedance  real  and 
imaginary  for  20%  noise. 
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Figure  61.  Cumulative  coherence  functions 
for  20%  noise. 
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